Phase synchronization and topological defects in inhomogeneous media 
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The influence of topological defects on phase synchronization and phase coherence in two- 
dimensional arrays of locally-coupled, nonidentical, chaotic oscillators is investigated. The motion 
of topological defects leads to a breakdown of phase synchronization in the vicinities of the defects; 
however, the system is much more phase coherent as long as the coupling between the oscillators 
is strong enough to prohibit the continuous dynamical creation and annihilation of defects. The 
generic occurrence of topological defects in two and higher dimensions implies that the concept of 
phase synchronization has to be modified for these systems. 
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The rich collective behavior, including mutual entrain- 
ment and self-synchronization, seen in systems of coupled 
oscillators has been a stimulus for the long-standing in- 
terest in these systems (see, for example, Refs. [Q, |[ |^ 
and references therein) . Recently, attention has turned to 
the study of coupled chaotic oscillators and, in particular, 
to the phenomenon of phase synchronization. Provided 
that the phase can be defined |, |], two coupled non- 
identical chaotic oscillators are said to be phase synchro- 
nized if their frequencies are locked 0, |. This appears 
to be a general phenomenon and it has been observed 
in such diverse systems as electrically coupled neurons 
[0, H , biomedical systems , chemical systems , and 
spatially extended ecological systems to name only 
a few. Moreover, the potential role of phase synchroniza- 
tion in brain functions has been recognized 

Much of the theoretical analyses of phase synchroniza- 
tion have been carried out on systems consisting of two 
locally coupled oscillators or many globally coupled 
oscillators [Q. Large one-dimensional chains of locally 
coupled chaotic oscillators have been investigated very re- 
cently ||l5, 16, 1^. Here, we address the question whether 
phase synchronization can persist in higher spatial di- 
mensions where topological defects can play a central 
role. We show that the existence of topological defects 
can lead to a breakdown of global phase synchronization 
in two-dimensional arrays of nonidentical chaotic oscilla- 
tors. While most of the medium may remain phase syn- 
chronized, oscillators close to moving topological defects 
have a different frequency. Despite this fact, the phase 
coherence of the system is higher than in systems with- 
out topological defects. The transition to phase synchro- 
nization via phase clustering observed in one-dimensional 
systems jl|, ^ is not found in our simulations on two- 
dimensional systems; instead, a transition involving point 
defects occurs. 

Point topological defects in two-dimensional homoge- 
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neous oscillatory media are associated with the appear- 
ance of spiral waves |2^, ^ . The phase field 0(r , t) of 
a medium with a spiral wave contains a point topological 
phase defect in the spiral core such that 



^^V0(r,t) -dl^int. 



(1) 



where nt is the topological charge of the defect A 
topological defect corresponds to a point in the medium 
where the local amplitude is zero and the phase is not 
defined. For periodic boundary conditions, the net topo- 
logical charge of the medium is zero. For identical chaotic 
oscillators pa, spatially coherent spiral dynamics can 
still exist implying a phase locking of the oscillations . 

To illustrate the phase synchronization properties in 
two-dimensional networks of nonidentical chaotic oscilla- 
tors, we consider an L x L array of locally coupled Rossler 
oscillators with periodic boundary conditions [E5[ , 



ax(r,t) 
dt 



R(x(r,i)) 



K ^ (x(f,t) 



x(r,i)), (2) 



where Ri = —Lo[r)x2 — ^3, R2 — LLi[r)xi + Q.2x2, R3 = 
X1X3 — 5.9x3 + 0.2. The sites of the lattice are labelled by 
r, K is the coupling constant and iV(r) is the set of the 
four next nearest neighbors of site r. The phase angles 



of the oscillators are given by 0(r, t) = arctan 



E2(r,t) 

Xi(T,t) 



To 



obtain the results presented below, we take L = 64 and 
choose the a;(r)'s randomly from a uniform distribution 
on the interval [0.95, 1.05], ensuring that the system is in 
the chaotic regime. 

First, we consider K = Kq = 0.1. For different initial 
conditions leading to stable patterns, basically two sce- 
narios can be observed. For homogeneous initial condi- 
tions or initial conditions with small inhomogeneities, no 
topological defects are created and the system evolves to 
a target pattern similar to the one in Fig. |l] (left panel). 
In the case of larger inhomogeneities in the initial con- 
ditions, a number of topological defects with = ±1 is 
created initially. Some pairs of topological defects with 
opposite charge are quickly annihilated until only a small 
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FIG. 1: Snapshots of the phase field, all for the same realiza- 
tion of uj{r). From left to right: K = 0.05 without topological 
defects, K — 0.05 with two surviving topological defects, and 
K — 0.0419. In the latter case, topological defects are gener- 
ated dynamically in pairs such that the number of topological 
defects fluctuates around a value of 8 pairs. 
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number is left generating a spiral pattern similar to that 
in Fig. |l| (middle panel). The topological defects are not 
necessarily stationary and for very long times it is possi- 
ble that all defects could disappear through further an- 
nihilation events; however, the motion of the surviving 
topological defects is very slow compared to the initial 
annihilation processes and no further annihilation events 
were observed on our long simulation time scale (10** - lO'"^ 
spiral revolutions). These two scenarios persist provided 
K is not too small. 

When there are no topological defects, the system 
is phase synchronized. The occurrence of phase syn- 
chronization for coupled Rossler oscillators is usually at- 
tributed to the high degree of phase coherence of the 
attractor of a single Rossler oscillator [|l| although this 
might not be a sufficient condition ||^. Phase coher- 
ence means that (/^ dT{T]{t)7]{t + T))t)i/2 ^ cjo, where 

d(j)/dt = Wo + '7(0; 4' is the unbounded phase and (• • ■)t 
signifies an average over t. If rj is (5-correlated Gaussian 
noise this condition reduces to ct ^ cjo and corresponds 
to a very narrow peak of width at ojq in the power spec- 
trum of xi{t). However, generally temporal correlations 
in ?7 exist. These correlations determine the speed of con- 
vergence of the time average ui{T) = T~^{<j){T) — 0(0)) 
towards wq. For the system of coupled Rossler oscillators 
considered here, the speed of convergence — as measured 
by the standard deviation of the ensemble distribution 
of w(T) — scales as 1/T as shown in Fig. ^. This is 
in contrast to what one would expect if the ry variables 
were independent where one obtains \/^/T scaling since 
uj{T) — ojQ = 1/r jj" dtrjit) and the standard deviation of 

dtri(t) scales with \/T. The observed variation with 
T is much faster, implying an extremely high degree of 
phase coherence, even in the coupled system. The origin 
of this effect resides in the shape of the local attractors. 
They are two-banded with distinct frequencies oji and uj2 
such that wq = (wi -|-a;2)/2. Moreover, each deviation 
from oji is followed by a deviation from ijj2 of similar am- 
plitude but with opposite sign, and vice versa. Hence, 



FIG. 2: Log-log plot of the standard deviation acj{T) of the 
cD(r) distribution in the phase synchronized state describing 
the speed of convergence towards wo- The solid line with slope 
-1 is plotted to guide the eye. Note that To — 27r/aJo depends 
slightly on K. 



the deviations of uj(T) from luq should be dominated by 
the switching between the bands for not too large val- 
ues of T. Thus, the maximal deviation from too is given 
by the term {uji — W2)/(2T), reproducing the observed 
scaling. The two-banded structure can also be found in 
the ensemble distribution of Cd{T). There, two distinct 
peaks can be seen for certain T values consistent with 
the fluctuations around the 1/T scaling in Fig. |[ 

If topological defects are present global phase synchro- 
nization is no longer guaranteed. Since a topological de- 
fect corresponds to a point in the system where a phase 
variable cannot be defined, this implies that the concept 
of phase synchronization can only be applied to parts of 
the system which exclude topological defects. For dis- 
crete lattices of oscillators considered here, this does not 
pose any problems because the defect is located almost 
surely between lattice sites. Thus, the presence of topo- 
logical defects does not exclude global phase synchroniza- 
tion a priori. We indeed observe phase synchronization 
for quasi-stationary topological defects which explains 
the occurrence of phase synchronization in the presence 
of topological defects described in Ref. . However, the 
motion of topological defects destroys phase synchroniza- 
tion locally. As shown in Fig. ^, the motion of a topo- 
logical defect leads to a distortion of nearby local orbits. 
This in turn influences their local frequencies such that 
the system is no longer phase synchronized as can be de- 
duced from the corresponding loq distribution in Fig. ^ 
pTf . In particular, the "outliners" in this distribution 
are located in the vicinity of the topological defect (s). 

For homogeneous oscillatory media it is well known 
that the time scales for the motions of topological de- 
fects depend on the system parameters, ranging from 
very slow motion to rapid dynamics 128] . Inhomogeneity 
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FIG. 3: Trajectory of x-i{t) (black) for a fixed oscillator and 
the distance Ar{t) (white) of this oscillator from the moving 
topological defect for K = 0.06. The two temporal behaviors 
are correlated, demonstrating the influence of the topological 
defect's motion on the local orbit. 
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FIG. 4; The uo histogram h{uJo) of the oscillators for 
K = 0.06 with two surviving topological defects in the sys- 
tem. One defect moves causing the breakup of global phase 
synchronization. The exact shape of the distribution depends 
on the w(r) realization. Negative and/or positive "ouliners" 
can exist and correspond to oscillators in the vicinity of the 
moving topological defect. 



provides another mechanism for defect motion and, in 
the present case of a network of non-identical Rosslcr os- 
cillators, this is the dominant mechanism. For a network 
of identical Rossler oscillators, the topological defects are 
quasi-stationary on the time scales considered here; thus, 
the observed breakup of global phase synchronization for 
this system is due to the quenched disorder in the net- 
work. 

Although global phase synchronization is lost in the 
presence of moving topological defects, the fluctua- 



tions of the instantaneous period in the medium are 
greatly reduced. While the phase coherence criterion, 

(/.^ '^'''^''(*)'?(* '^))ty^^ ^0, is satisfied to a high 
degree of accuracy as in the case without topological de- 
fects, the average amplitude of the fluctuations in the 
instantaneous periods T{r,n) = 27r/w(r, n) as measured 
by5=((5(r))r)^ with 



Sir) 



(r(r,7i))„ 



v/((T(r,n)-(r(r,n))„)2 



(3) 



are significantly different. Here, T(r, n) is the time 
needed for the nth rotation at site r and (• • is the 
average over all realizations of cj{r). Figure]^ (left panel) 
shows that the value of S is much larger above a certain 
coupling Kc when (more) topological defects are present 
in the system, corresponding to a higher degree of coher- 
ence. An analysis of S{r) shows that high values exist, 
especially in the vicinity of quasi-stationary topological 
defects (see right panel of Fig. ||). This is due to the 
influence of topological defects on the shapes of the local 
orbits: they are close to periodic limit cycles. This is in 
accord with observations in for identical oscillators 
although the "effective" spatial period doubling cascade 
found there is absent in the present case. 

The system maintains memory of its initial condition 
above Kc ~ 0.043 leading to the different values of S. For 
single a;(r) realizations, different scenarios, with either 
no defects or with defects present, can be observed for 
A'o > K > Kc For K below Kc, however, a qualita- 
tively new behavior occurs. Topological defects are now 
generated and annihilated continuously by the dynamics 
which is responsible for the loss of memory of the initial 
condition. The number of topological defects fluctuates 
around a mean value which increases with decreasing K . 
This is expected because K should control the character- 
istic length scale in the system. Moreover, the topologi- 
cal defects move on average much faster than above Kc, 



similar to defect mediated turbulence |31 . The creation 



and annihilation of topological defects as well as the fast 
motion lead to strong distortion of the local orbits af- 
fected by the moving defects . This makes it impossible 
to define a proper phase variable for most of the oscilla- 
tors in the network. Consequently, the concept of phase 
synchronization can no longer be applied. 

In Refs. ||l5|, the transition to phase synchroniza- 
tion via the route of phase clustering and merging of 
clusters was investigated for chains of oscillators. Liu, 
et al. concluded that phase clustering should be 

more prevalent than full phase synchronization, espe- 
cially in networks of coupled neurons. However, the two- 
dimensional system considered here shows a different be- 
havior for increasing K starting from a value below Kc- 
We observe a transition from a state where topological 
defects are generated continuously by the dynamics to a 
partial phase synchronized state with moving topological 
defects, similar to the case with surviving point defects. 
This partial phase synchronized state is characterized by 
a large number of phase synchronized oscillators and a 
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FIG. 5: Left: 5 as a function of K for homogeneous (star) 
and inhomogeneous (diamonds) initial conditions. The results 
were obtained from an average over five realizations of the 
random uj{r) field. The phase coherence is higher when (more) 
topological defects are present in the system for K > Kc ~ 
0.043. Below Kc point defects are created and annihilated 
dynamically. Right: S(r) for the same u>{r) realization as in 
Figs. ^ and|^. Dark areas correspond to high values of ^(r). 
The locations of the defects are given by the two arrows. The 
lower left defect is the moving defect. 



small number of oscillators close to the point defects with 
mutually different frequencies. Thus, only one cluster 
with more than one oscillator exists in contrast to the 
chain geometry. The existence of topological defects in 
two dimensions and higher dimensions suggests that the 
transition described here should be prevalent for general 
networks of chaotic oscillators. 

We have shown that the occurrence of topological de- 
fects in two-dimensional arrays of nonidentical chaotic 
oscillators has a two-fold effect: moving topological de- 
fects lead to a breakup of global phase synchronization, 
while all topological defects — especially quasi-stationary 
ones — increase the phase coherence of the system signifi- 
cantly. Since moving and quasi-stationary topological de- 
fects are present in many circumstances, the apparently 
contradictory scenario of increased phase coherence but 
loss of global phase synchronization occurs. 
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